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I. INTRODUCTION 



The phenomenology of neutralino dark matter has been studied extensively in the Mini- 
mal Supersymmetric extension of the Standard Model (MSSM) HJ. This model incorporates 
the same gauge group as the Standard Model and the supersymmetric extension of its parti- 
cle content. The Higgs sector is slightly modified as compared to that of the Standard Model: 
the MSSM requires two Higgs doublets Hi and H 2 in order to give mass both to down- and 
up-type quarks and to cancel anomalies. After Electro- Weak Symmetry Breaking (EWSB), 
the physical Higgs fields consist of two charged particles and three neutral ones: two scalar 
fields (h and H) and one pseudoscalar (A). The Higgs sector is specified at the tree level by 
two independent parameters: the mass of one of the physical Higgs fields and the ratio of 
the two vacuum expectation values, usually defined as tan/3 = v 2 /v x =< H 2 > / < Hi >. 
The supersymmetric sector of the model introduces some other free parameters: the mass 
parameters M x , M 2 and M3 for the supersymmetric partners of gauge fields (gauginos), 
the Higgs-mixing parameter \i and, in general, all the masses of the scalar partners of the 
fermions (sfermions). 

In the MSSM it is generally assumed that the gaugino masses are equal at the grand 
unification scale Mqut'- M^Mqut) = m i/2 and hence are related at lower scales by 

Mi : M 2 : M 3 = a x ■ a 2 : a 3 (1) 

where the «j (i=l,2,3) are the coupling constants of the three Standard Model gauge groups. 
The neutralinos are mass-eigenstate linear superpositions of the two neutral gauginos (7 and 
Z) and the two neutral higgsinos (Hi and H 2 ) 

X = 017 + a 2 Z + a 3 Hi + a 4 H 2 . (2) 

The neutralino sector depends, at the tree-level, on the following (low-energy) parameters: 
Mi = (5/3) tan 2 6y/M 2) M 2 ~ 0.8 m^, fi and tan (5. Neutralino properties are naturally 
discussed in the (m^, fi) plane, for a fixed value of tan/3. As an example, in Fig.l the lines 
of constant mass for the lightest neutralino (m x ) and constant gaugino fractional weight 
(P = af + a 2 ,) are plotted in the (mi/ 2 , fi) plane for tan/3 = 8. We observe that the mass of 
the lightest neutralino increases from the bottom left to the top right, while the neutralino 
composition changes from higgsino dominance in the top-left region of the plane to gaugino 
dominance in the bottom-right. The regions forbidden by accelerator data are also displayed 
in Fig.l. 

The low-energy MSSM scheme is a purely phenomenological approach, whose basic idea 
is to impose as few model-dependent restrictions as possible. In this approach the lightest 
neutralino is a favourite candidate for cold dark matter. This scheme has been employed ex- 
tensively in the analysis of the size and the relevance of various possible signals of neutralino 
dark matter: direct detection , signals due to neutralino annihilation in celestial bodies, 
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namely the Earth and the Sun |5]||, and signals from neutralino annihilation in the galactic 
halo [0. The MSSM provides a useful framework in which neutralino phenomenology may 
be analysed without strong theoretical prejudices which could, a posteriori, turn out to be 
incorrect. This scheme is also frequently employed in analyses of the discovery potential of 
future accelerators ||. 

At a more fundamental level, it is natural to implement this phenomenological scheme 
within the supergravity framework fl9|-|TT|] . One attractive feature of the ensuing model is the 
connection between soft supersymmetry breaking and EWSB, which would then be induced 
radiatively. The essential elements of the model are described by a Yang-Mills Lagrangian, 
the superpotential, which contains all the Yukawa interactions between the standard and 
super symmetric fields, and by the soft-breaking Lagrangian, which models the breaking of 
supersymmetry. Here we only recall the soft supersymmetry breaking terms 

- C S oft= ^m-|0i| 2 

i 

+ { [A l ab h l ab L a HxR b + A^QaH^ + A u ab h u ab Q a H 2 U b + h.c] - B^H X H 2 + h.c.} 
+ ^2 Mi(XiXi + XiXi) (3) 



where the & are the scalar fields, the Aj are the gaugino fields, Hi and H 2 are the two Higgs 
fields, Q and L are the doublet squark and slepton fields, respectively, and U, D and R 
denote the S"?7(2)-singlet fields for the up-squarks, down-squarks and sleptons. In Eq.(|), 
rrti and Mj are the mass parameters of the scalar and gaugino fields, respectively, and A 
and B denote trilinear and bilinear supersymmetry breaking parameters, respectively. The 
Yukawa interactions are described by the parameters h, which are related to the masses of 
the standard fermions by the usual expressions, e.g., m t = h l v 2 . 

The supergravity framework is usually implemented with a number of restrictive assump- 
tions about unification at Mqtjt'- 

i) Unification of the gaugino masses: M^Mgut) = ^1/2, 

ii) Universality of the scalar masses with a common mass denoted by m : iti^Mgut) 

= m , 

iii) Universality of the trilinear scalar couplings: A\Mgut) = A d (McuT) = A u (Mgut) 

= A m . 



These conditions have strong consequences for low-energy supersymmetry phenomenology, 
and in particular for the properties of the neutralino as dark matter particle. Typically, the 
lightest neutralino is constrained to regions of gaugino dominance, that entail a large relic 
abundance (in wide regions of the parameter space fi x h 2 exceeds the cosmological upper 
bound) and a small direct detection rate for neutralino dark matter. Indirect signals from 
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the neutralino, such as high-energy neutrinos from the Earth and Sun, and the products of 
annihilation in the halo, are practically undetectable [11]. 



The above assumptions, particularly ii) and iii), are not very solid, since universality 
may occur at a scale higher than Mqut, i-e-, the Planck scale or string scale [|1^], in which 
case renormalization above Mqut may weak universality in the m;, e.g., between scalars in 
5 and 10 representations of SU(5) |fT3| . Moreover, in many string models the m^'s are not 
universal even at the string scale. 



In a number of recent works ]T4|,[T5|, deviations from some of the unification conditions 
have been considered. In particular, in Ref. [Q phenomenological consequences for neu- 
tralinos of a relaxation of assumption ii) have been analysed in the regime of large values 
of tan p. It has been shown that deviations from condition ii) may entail a changeover in 
neutralino composition from a gaugino-like state to a higgsino-like state (or at least to a 
higgsino-gaugino mixed state), with important consequences for neutralino phenomenology. 



In this paper, we first explore, over the full range of tan (3, the various scenarios which 
may occur when condition ii) is relaxed, with an approach which is similar to the one adopted 
in the large-tan (3 analysis of Ref. fll4j[ . We then discuss in detail the ensuing consequences 
for neutralino dark matter, with particular emphasis for its direct detection. 



In the following, we first discuss which constraints can be applied to the parameters 
when specific physical requirements are imposed. In Sect. II, we summarize the conditions 
implied by radiative EWSB and define the type of departure from universality examined in 
this paper. Then, in Sect. Ill we establish some upper bounds on the supergravity param- 
eters by requiring that radiative EWSB does not occur with excessive fine tuning among 
different terms. In Sect. IV we analyse in detail the constraints due to the requirement that 
EWSB takes place radiatively. Subsequently, in Sect.V cosmological constraints, derived 
from the evaluation of the neutralino relic abundance, are discussed. Other constraints, 
from experimental data on b — > sj processes and on the mass of the bottom quark w,b, are 
applied in Sect. VI. In Sect. VII the effects of these various constraints are first displayed in 
the (mi/2, m ) plane for fixed tan/3 and A Q , and then shown in the (m.1/2, yu) plane, which 
provides the most useful representation for discussing neutralino phenomenology. We recall 
some specific properties of the neutral Higgs bosons in Sect. VIII. Finally, in Sect. IX event 
rates for direct detection of neutralino dark matter are discussed. Conclusions are presented 
in the last Section. 



II. RADIATIVE EWSB 

We recall that the tree-level scalar potential for the neutral Higgs fields may be written 
in the form 

V = (M 2 Hi + /i 2 )|#i| 2 + {M 2 H2 + fi 2 )\H 2 \ 2 - Bn{HxH 2 + h.c.) + quartic D terms. (4) 
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The parameters of this potential must obey the following physical conditions: 

Sin2/3= M Sl+ lf + 2^ < 5 > 

Ml = ^k-M^ff _ 

z tan 2 (3 - 1 r w 

Afl = M ^ + ^4 2 + 2/i 2 > . (7) 

Here M4 is the mass of the CP-odd neutral Higgs boson (see Sect. VIII below), and eq.(f^) 
must in fact be strengthened to M4 > {Mj^ibi where (Ma)w is the experimental lower 
bound fl6|]. For instance, for tan/3 > 3, (Ma)/& ~ 45 GeV. Notice that the sign of ji is 



defined according to the convention of reference |]]. We remark that although Eqs.([|-[FD are 
expressed at the tree level, in our actual calculations 1-loop corrections to V [|17] have been 
included. The M#.'s (as well as the sfermion and the gaugino masses and the parameters A, 
B and //) evolve from the Mqut scale down to the Mz scale according to the Renormalization 
Group Equations (RGE's). This is how Eq.(B) may be satisfied, even if M Hl and M H2 are 
equal at M GUT . 

In this work we consider deviations from universality in the scalar masses at Mqut-, 
which split Mh 1 from Mh 2 - This effect is parameterized as 

M 2 (M GUT )=m 2 (l + 5 i ) . (8) 



The parameters 5,i which quantify the departure from universality for the M\, will be varied 
in the range (— 1,+1), but are taken to be independent of the supersymmetry parameters. 
This is an Ansatz, since, when evolving the scalar masses from the unification scale (Planck 
scale or string scale) to the GUT scale Mqut, the deviation parameters are in general 
functions of all the supersymmetry parameters fl8|| . 



Following a common procedure, Eq.fl^) is used to replace the parameter B by tan/3. Thus 
the set of independent parameters becomes mi/2, itiq, Aq, tan/3, and fi 2 is given in terms of 
these parameters by Eq.(^), suitably corrected by 1-loop effects: only the sign of [i remains 
undetermined. Obviously, values of /z 2 are accepted only if they exceed the experimental 
lower bound fif b , which is derived from the lower limit on the chargino mass ||16|| : — 45 
GeV. 

We have solved the RGE's using the 1-loop beta functions including the whole super- 
symmetric particle spectrum from the GUT scale down to Mz, neglecting the possible effects 
of intermediate thresholds. Two-loop and threshold effects on the running of the gauge and 



Yukawa couplings are known not to exceed 10% of the final result [19|. While this is of 



crucial importance as far as gauge coupling unification is concerned |T3 , it is a second-order 
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effect on the evolution of the soft masses. Since neutralino properties are studied over a wide 
range of variation for the high scale parameters, such a degree of refinement is not required 
here. 



In order to specify the supersymmetry phenomenology, boundary conditions for the gauge 
and Yukawa couplings have to be specified. Low-scale values for the gauge couplings and 
for the top-quark and the tau-lepton Yukawa couplings are fixed using present experimental 
results. In particular, we assign for the top mass the value m t = 178 GeV pD| . In addition, 
we require the unification of the bottom and tau Yukawa couplings at the GUT scale, as 



would be suggested by a unifying group that includes an S77(5)-like structure pi 



The values of Mh x and Mh 2 at the Mz scale, obtained from the RGE's, may be param- 
eterized in the following way: 

Mff. = aim\ /2 + biUil + CiA%m% + diA m m 1/2 ■ (9) 



(Notice that, in our notation, all running quantities written without any further specification 
are meant to denote their values at Mz-) The coefficients in the expression @ are functions 
of tan (3 and of the <Vs. They are displayed in Fig. 2 (a,b) for the case of universal scalar 
masses, {i.e., 5i = 0). The coefficients for Mfj 2 turn out to be very stable as functions of 
tan/3, except for small tan/3. More precisely, a 2 ~ —2.5 for tan/5 > 4 with all the other 
coefficients much smaller (of order 0.1). As far as Mfj is concerned, whereas C\ and d\ are 
again very stable (of order 0.1), a\ and bi vary rapidly as functions of tan/3. This property 
of ai and b\ is due to the very fast increase of h b for increasing tan f3. 

When a departure from mo universality is introduced, the coefficients in Eq.(H), except 
for cii and a 2 , become functions of the parameters 5f. b±, c\ and d\ depend on 5± and b 2 , c 2 
and d 2 on 5 2 . Whereas the fej's are rapidly-increasing functions of the <5j's, the Cj's and the 
di§ are rather insensitive to these parameters. 

Stringent constraints on the parameters mi/ 2 ,mo, Aq and tan/3 follow from the request 
that the M^.'s, evaluated from Eq.(^), satisfy Eqs.(|G|-[7p. Explicitly, we require that /i 2 and 
M\, given by the expressions 

/i 2 = - — 2^ — - «2 tan 2 (3)m\, 2 + (bi - b 2 tan 2 (3)m 2 Q + 

M 2 

(ci - c 2 tan 2 (5)A\m^ + {dx - d 2 tan 2 l3)AQm mi/ 2 } — 

M 2 

= Jim 2 /2 + J 2 ml + J s Alml + J 4 Aom m 1/2 ^ (10) 

M\ = (ax + a 2 + U x )m\ l2 + (bt + b 2 + 2J 2 )m 2 + 

(ci + c 2 + 2J 3 )A 2 m 2 + (dt + d 2 + 2J i )A m m 1/2 - M| 
= K x m\ l2 + K 2 m\ + K 3 A 2 m 2 + K 4 A m m l/2 - M% (11) 
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satisfy the conditions: fi 2 > /i 2 b , Ma > (Ma)z& mentioned earlier. 

The coefficients Jj and Ki in Eqs. (|10| , [TTD are plotted as functions of tan/3 in Fig.2 (c,d) 



for the case of m universality. In Fig. 2c we notice that all the Jj's are positive, with Ji 
dominating over the others: for tan /3 > 4, one has J% ~ 2.4. As far as the coefficients K$ are 
concerned, we see in Fig. 2d that only two of them are sizeable: K\ and Ki- They are both 
decreasing functions of tan/3, with K\ > Ki- At very large tan/3 these coefficients become 
very small, and K 2 even becomes negative (but still small in magnitude) at tan (3 > 50. 

In the case of non-universality, the coefficients Jj and Ki, except for Ji and Ki, become 
functions of the parameters Si. We will see in Sect. IV that many important features of the 
supersymmetry parameter space depend on the signs of the two coefficients J 2 and K 2 . We 
show in Figs. 3 and 4 how their signs depend on the values of the <5j's. In Fig. 3 the lines 
J 2 = are plotted in the (5 2 , <5i) plane for a few values of tan/3: for each value of tan/3, 
J 2 is negative in the region above the relevant J2 = line and positive below. Similarly, in 
Fig. 4 the K 2 = lines are displayed in the same (5 2 , 5i) plane at fixed tan/3: K 2 is negative 
above the K2 = lines, and positive below. 



We now make a few comments related to Eq. ([TT|) , since the value of Ma plays a very 
crucial role in a number of important neutralino properties. This is due to the fact that 
many physical processes involving neutralinos are mediated by the neutral Higgs bosons. 
Thus the value of Ma determines the size of the relevant cross sections both through Mi- 
dependence in propagators and, in an implicit way, through the couplings of the h and H 
bosons to quarks and to the lightest neutralino x ( see Sect. VIII). As a consequence, a small 
value of Ma has the effect of enhancing the magnitude of the relevant cross sections. 

What values of Ma do we obtain from Eq. (|TTD ? Because of the properties of the coeffi- 
cients Ki previously analysed, Ma turns out to be a rapidly-decreasing function of tan (3. In 
Fig. 5, Ma is displayed at the representative point m = 50 GeV, mi/ 2 = 200 GeV (1-loop 
corrections to Ma have been included in the calculation). One notices that Ma is OiMz) 
for tan/3 > 45. This feature provides one of the most appealing scenarios for neutralino 
phenomenology. 



III. CONSTRAINTS DUE TO THE ABSENCE OF FINE TUNING 



Before we exploit fully the two constraints /i 2 > fif b , Ma > (M^b to restrict the pa- 
rameter space, we apply the general criterion that the expression ( |10D is satisfied without 
excessive tuning among the various terms |22]JT0|] . In radiative EWSB the physical value of 
Mz, which sets the EW scale, may be written as 



Ml 



2(Jim 2 /2 + J 2 ml + J 3 Alml + J 4 A m m 1/2 - /J 2 



(12) 
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Accidental compensation (fine tuning) among different terms in Eq.(|l"2"D may occur. We 
explicitly require the absence of too-strong fine tuning, i.e., cancellations among exceedingly 
large values of the parameters m^, m , A and fi. Denoting by rjf a parameter which 



quantifies the degree of fine tuning, we require p2| that 



AMf 



< rj f 



xf 



(13) 



where Xi denotes any of the previous parameters. For instance, for A = 0, Eq.(|13"|) provides 
the following conditions 

m '/ 2< 2^ M| ' m ° < 2j!^ M| ' ^ < " (64 ° GeV)2 (M) 



where in the last approximate equality we have taken 77/ = 100, which means that we allow 
an accidental compensation at the 1% level. The upper bound on m depends on tan/3 and 
the SiS, whereas that on varies only with tan/3 (because of the nature of the Ansatz 
(§): see the comment after Eq.(|8|)). 

For the sake of illustration, we give some numerical examples, taking again rjf = 100. 
For tan (3 = 8, we have, for S\ = S 2 = 0, mi/2 < 400 GeV, mo < 1.5 TeV. For two other pairs 
of values of the S^s, which will be discussed later on, we obtain mo < 2.4 TeV for 5\ = —0.2, 
5 2 = 0.4 and m < 3.0 TeV for Si = -0.8, 5 2 = 0.2. At tan/3 = 53 we have m l/2 < 415 GeV 
and m < (1.7 — 1.9) TeV, depending on the values for the S^s. These inequalities imply for 
the neutralino mass m v < 170 GeV. 

In the following, when graphical representations for the parameter space are shown, 
we display no-fine-tuning upper bounds obtained from the general expression ([13]) with 
rjf = 100. These upper bounds are denoted by dashed lines in Figs. 9-14. 



IV. CONSTRAINTS DUE TO RADIATIVE EWSB 



The EWSB constraints are given by the set of Eqs.(|-[7]), or equivalently by Eqs. (|l0H Tl|), 
together with the conditions /i 2 > jif b and Ma > (Ma)u>- From these equations the values 
of m and (or ji and m^) are constrained and thus some domains in the (m 1( / 2 , m ) 
or (7711/2, aO planes can be excluded. Let us start this discussion by analyzing the condition 
Ma > {Ma) ib, with Ma given by Eq. (|TT|) ■ For the sake of simplicity, we put A = for the 
moment. To discuss the role of Ma > {Ma)h> in placing bounds on mi/2 and m , we first 
rewrite it explicitly as 

K x m\, 2 + K 2 m\ > M% + {M A f lb . (15) 
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The nature of this quadratic form in the (my 2 , m ) plane obviously depends on the signs 
of the two coefficients K\ and K 2 . As we have seen in Sect. II, it turns out that, whereas 
Ki is always positive, the sign of K 2 depends on the values of tan/3 and of the S^s. Two 
different situations may occur, depending on the sign of K 2 . In the case K 2 > the region 
allowed by (|T5| ) is the one above an elliptical branch centered in the origin of the (mi/2, ttiq) 
plane. Therefore, both parameters vn\j 2 and mo are bounded from below. When K 2 < 0, 
the region allowed by Eq.flToD is the one between the rrii/ 2 axis and an upward-moving 
hyperbolic branch. Thus, whereas mi/ 2 is still bounded from below, mo is now constrained 
from above. The upper bound on m is particularly stringent when K 2 is large and negative 
and Ki is not large. This occurs, for instance, at very large values of tan/3 in the case of 
mo universality. 

This discussion may be extended straightforwardly to the case Aq 7^ 0. In this case the 
constraint M4 > (M^)^ may be written explicitly as 

K im \ /2 + K 2 ml + K 3 A 2 m 2 + K 4 A m m l/2 > M% + (M A f lh . (16) 



The nature of this quadratic form depends on the sign of its determinant. When this 
determinant is positive, an elliptical branch in the (m^, mo) plane provides lower bounds 
on the two variables. On the other hand, a negative determinant entails an upward-moving 
hyperbolic branch which places an upper bound on m . These branches are part of conies 
whose axes are somewhat tilted with respect to the (m^, m ) axes. 

Similar implications follow from the constraint fi 2 > fi 2 b , which may be written explicitly 
as (for A = 0) 

M 2 

J im ? /2 + J 2 mg>-^ + /4. (17) 



This quadratic form may be discussed in much the same way as the one in Eq.(|T^). From 
the properties seen in Sect. II it turns out that the coefficient J\ is always positive, whereas 
the coefficient J 2 is positive in the universal case, but may be negative when deviations 
from m universality are introduced. Thus it follows that the condition fi 2 > fif b puts lower 
bounds on mi/ 2 and either lower or upper bounds on m , depending on the sign of J 2 (due 
to analytic properties identical to those discussed previously below Eq.(|l|)). The condition 
> se ^ s a very stringent upper bound on m , whenever J 2 is negative and large in 
magnitude. The extension to the case A ^ may be repeated here in a way similar to the 
above discussion for Eq. (|16|) . 

Thus we have seen that two important constraints, ji 2 > fif b and Ma > (M^jib, are at 
work in bounding and mo, when EWSB is required to occur radiatively. When J 2 and 
K 2 are positive, the two conditions place lower bounds on and m . Similar constraints 
are established by the requirements that also the sfermion masses and m x satisfy the relevant 
experimental bounds. These last conditions are not explicitly discussed here, but they are 
taken into account in our evaluations. 
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It is worth emphasizing that the most dramatic impact of the conditions fi 2 > fif b and 
Ma > (Ma)w over the parameter space occurs when either J 2 or K 2 (or both of them) are 
negative. Under these circumstances, as we have seen above, /i 2 > fif b and Ma > (Ma)w may 
place stringent upper limits on m , bounding the neutralino parameter space considerably. 
Which of the two conditions prevails over the other depends on the specific regions of the full 
parameter space and on the values of the <Vs. In Sect. VII we will illustrate the implications 
of these constraints in a few specific examples. 



V. COSMOLOGICAL CONSTRAINT 



Let us turn now to the evaluation of the neutralino relic abundance Q x h 2 and to the 
requirement that the lightest neutralino is not overproduced, i.e., Q x h 2 < 1. 

The neutralino relic abundance Q x h 2 is evaluated following the standard procedure 
|2"B|^Z6], according to which VL x h 2 is essentially given by Vt x h 2 oc< a^v >^J, where 
< &aiu\V >int is the thermally-averaged annihilation cross section, integrated from the freeze- 
out temperature to the present temperature. The standard expansion < cr ann t> >= a + bx + ... 
may be employed, with x = T/m x , except at s-channel resonances (Z, A, H, h), where a more 
precise treatment has to be used for the thermal average ||24|| . In the evaluation of < cr ann t> > 



the full set of annihilation final states (// pairs, gauge-boson pairs, Higgs-boson pairs and 
Higgs-gauge boson pairs), as well as the complete set of Born diagrams are taken into ac- 
count [F2BJ. We recall that one of the largest contributions to the annihilation cross section is 
provided by diagrams with the exchange of the pseudoscalar Higgs boson A. (More relevant 
properties of the Higgs bosons are discussed in Section VIII.) We note that the constraint 
^x^ 2 ^ 1 is very effective for small and intermediate values of tan (3, but is not restrictive for 
large values of tan (3. The strong restriction in the former case comes from the large value 
of Ma implied by small and intermediate values of tan /3 (see Fig. 5) (also the couplings of 
A to x an d fermions are small for these values of tan/?). 

We show in Figs. 6-8 a few examples where Q x h 2 is given as a function of m x in the 
form of scatter plots. These scatter plots have been obtained by varying the parameters 
mo and m^j on a equally-spaced linear grid over the ranges 10 GeV < mo < 2 TeV, 
45 GeV < mi/2 < 500 GeV. Furthermore, we remark that all evaluations presented in this 
paper are for positive values of fi, since negative values of \i are disfavoured by the constraints 
due to nib and b —>■ sj processes (see Sect. VI). The configurations shown in Figs. 6-8 satisfy 
the constraints due to radiative EWSB, discussed previously. 

In Fig. 6 is shown the case tan/? = 8 and Si = 0. Here, as expected because of the 
intermediate value of tan/3, many neutralino configurations provide Q x h 2 > 1, whilst only 
a few give fi x h 2 < 1. (Also, Ma is large here because of sizeable values of K 2 (see Fig. 2d), 
which helps increase fl x h 2 .) An exception occurs when m x ~ Mz/2, since in this case the 
annihilation cross section is greatly enhanced due to the Z-pole contribution. 
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In Fig. 7 we display £l x h 2 in a case of non-universality (8\ = —0.2, 5 2 = 0.4, for definite- 
ness). It is easier to find Q x h 2 < 1 in this case, since here the departure from m universality 
implies a changeover of the neutralino composition from the gaugino dominance of the pre- 
vious example to higgsino dominance (this point will be elucidated in Sect. VII). This implies 
a larger x~X annihilation cross section and consequently a smaller relic abundance. Thus 
only a few neutralino configurations are excluded by the Q x h 2 < 1 condition. 

An example for Q x h 2 in the case of large tan f3 and 5% = is shown in Fig.8. We see that 
^x^ 2 — 1 i m P oses no constraint since, for this very large value of tan f3, annihilation cross 
sections are very large. 



VI. CONSTRAINTS FROM b -> S7 AND m b 



In the evaluation of the b — > 57 decay rate we have included the supersymmetric con- 
tributions arising from the charged Higgs loops and chargino loops given in Ref. P7|] . The 
Higgs term always adds to the Standard Model value and usually entails too large a value for 
the rate. On the other hand, the chargino contribution gives rise to a destructive interference 
for fi > (in our convention for the sign of //). At large tan/3 supersymmetric contributions 
may be sizeable: unless the destructive interference protects the decay rate, it can very 
easily be driven out of the present experimental bounds. In the light of this property, the 
positive jj, scenario appears to be the favourite one and, as already remarked, in this paper 
we only show results for this case. In comparing our predictions with observations we have 
taken into account that, as discussed in Ref. ||28|| , large theoretical uncertainties are present, 



mainly due to QCD effects. In particular, predictions depend very strongly on the choice 
of the renormalization scale, leading to an inaccuracy of order 25%. To account for this 
effect we have relaxed the experimental bounds of Ref. p9f by the same amount, keeping 
the renormalization scale fixed at the representative value of 5 GeV. Thus, our requirement 
is that the rate of b — > s'j decay falls into the range 0.8 x 10~ 4 < BR(6 — > 57) < 5.3 x 10 -4 . 

The supersymmetric corrections to the bottom mass include contributions from bottom- 



squark-gluino loops and from top-squark-chargino loops p0 |. In the present analysis, the 
bottom mass is computed as a function of the other parameters and required to be compatible 
with the present experimental bounds. Theoretical uncertainties in the evaluations of m& 
arise both from the running of the RGE's and from assumptions about Yukawa unification. 
Since our choice is to solve RGE's at the 1-loop level and without thresholds, we estimate 
an uncertainty of the order of 10% in our prediction for m^. In addition, a relatively 
small departure (see Ref. |Jl[]) from bottom-r Yukawa unification at the GUT scale may 
significantly change the bottom mass result. To take into account such uncertainties we have 
chosen to weaken the bounds on m b given in [^| by an amount of 10%. Thus we require m b 
to fall into the range 2.7 GeV < m b (M z ) < 3.4 GeV. 
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VII. ALLOWED REGIONS IN NEUTRALINO PARAMETER SPACE 



We discuss now in a few examples how the various constraints analysed in the previous 
Sections complement each other in shaping the allowed regions in the parameter space. We 
start with the (mu2, mo) representation, and later display our results in the (m.1/2, y) plane 
which provides the most useful representation for neutralino phenomenology. 



Let us first clarify a few graphical conventions adopted in our (mi/2, m ) and (m 1 / 2 , //) 
plots. Regions are left empty when at least one of the following constraints is not satisfied: 



i) experimental bounds on Higgs, neutralino and sfermion masses [|lq , [33| , ii) the x is the 
Lightest Supersymmetric Particle (LSP), iii) radiative EWSB and fi 2 > fif b , Ma > (Ma)/&. 
Regions forbidden by the cosmological constraint (Q x h 2 < 1) are explicitly denoted by dots 
and those disallowed by the b — > 57, nib constraints (but not by the previous ones) are 
denoted by crosses (crosses are displayed only in the (mi/ 2 , m ) plane, but not in the (m 1 / 2 , 
ti) plane, to simplify these plots). The allowed domains are denoted by squares when they 
satisfy Q x h 2 > 0.01, or by diamonds otherwise in the (mi/2, m ) plots. They are denoted by 
squares in the (mi/2, /i) plots, independently of the fl x h 2 value. To simplify the discussion, 
we first take A = 0. We comment on the A 7^ case at the end of this Section. 



As a first example, let us consider the representative point tan/3 = 8. For this intermedi- 
ate value of tan f3, the cosmological constraint is expected to be very effective in view of the 
arguments discussed in Sect.V. This is actually the case for universal m , when both K 2 and 
J 2 are positive (see Fig. 2), so that the conditions of radiative EWSB do not set any upper 
limit on mo (Fig.9a). The empty region in the lower part of these figures is forbidden by 
the experimental bound on m x . As shown in this figure, in wide regions (denoted by dots) 
Q x h 2 > 1. Thus the cosmological constraint places a very stringent upper bound on mo for 
"^1/2 > 150 GeV. However, for smaller values of mi/2, an allowed horizontal region extends 
up to m — 2 TeV. In fact, along this strip, m x ~ Mz/2 and then fl x h 2 < 1 is satisfied (see 
the discussion in Sect.V). 

Moving away from the universal point towards a region where J 2 is negative, we expect 
^ 2 > 1*% to be effective in placing a stringent upper bound on m . This is actually the case 
in the example shown in Fig. 10a, which refers to the representative point 5± = —0.2, 5 2 = 0.4 
(J 2 = —0.07). Here it is the bound fi 2 > fif b which provides the most stringent constraint 
in disallowing the large (empty) domain on the right side. Nevertheless, Q x h 2 < 1 is still 
effective in excluding an internal region that would otherwise be allowed (see the discussion 
below) . 

Keeping tan/5 = 8, we complete our discussion by considering the representative point 
5\ = —0.8,52 = 0.2 shown in Fig. 11a, which gives an example where J 2 is very small. The 
peculiarity of this example will become clear when we discuss the relevant situation in the 
(mi/2, /-0 plane, to which we now turn. 

The shape and general properties of the physical region in the (m.1/2, yu) plane are dictated 
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by the constraints previously derived, and they are determined most notably by J 2 . It is 
convenient to distinguish the two cases i) J 2 > and ii) J 2 < 0. For case (i) at fixed m^, 
H increases for increasing m with the consequence that the allowed physical region extends 
to the right of the m = line in the (^1/2, /•*) plane, allowing for the neutralino only 
a gaugino-dominated region. In the case (ii) (J 2 < 0), starting from the = line and 
increasing mo at fixed m.1/2, one moves to the left and then one may reach regions of sizeable 
higgsino-gaugino mixing or even of higgsino dominance. Case i) applies in particular to the 
case of m universality (Si = 0) for any value of tan/5. This is clear from Fig.3, which shows 
that in the (#2, Si) plane the origin is below any J 2 = line. An example of this situation 
is displayed in Fig. 9b (for tan/3 = 8). 

However, as we have seen in Sect. II, when the assumption of mo universality is relaxed, 
then J2, which in the universal case is positive and small, may very easily become negative 
and sizeable. In this case a changeover in neutralino composition from an originally gaugino- 
like state into a higgsino-like one occurs. This remarkable property, discussed in Ref. |TJj 



for large tan /3, is in fact valid over the whole range of tan f3, if the degree of non-universality 
is increased for decreasing tan/3. An example of case ii) (J2 < 0) is shown in Fig. 10b, where 
the allowed region extends widely into the higgsino region. It is instructive to compare Fig. 9 
with Fig. 10. Looking at sections a) of these figures, we notice that changing the values 
of the S^s from the set Si = to the set Si = —0.2, S 2 — 0.4 relaxes substantially the 
cosmological constraint. Parts b) of these figures provide the explanations for this feature. 
In fact, whereas in the former case the neutralino is mainly a gaugino, in the latter case 
X is higgsino-like or mixed. As we already remarked, this implies an increase of the x~X 
annihilation cross section and a reduction of the relic abundance. The physical region also 
displays an extension to the right, in the example of Fig. lib, but here the effect is very tiny, 
due to a very small J 2 and to the severe upper bound on m for 7711/2 > 180 GeV. This is 
the first case to show a very marked (mi/2, \i) correlation. 

Now we turn to the case of large tan/3, where new features appear. First, the Ma > 
(Ma) ib condition is no longer protected by large values of Ki, and may become effective in 
restricting the parameter space. Secondly, the m& and b — * 57 conditions are now rather 
stringent over large domains and not only occasionally relevant as in the smaller tan/3 
cases. Thirdly, the cosmological constraint is usually overwhelmed by the other conditions. 
In Figs. 12a, 13a, 14a we have, for tan/3 = 53, the following sequence of examples, i) 
Si = 0, 82 — (Fig. 12a): here K 2 < 0, J 2 > 0, and since K 2 is negative and sizeable in 
magnitude, the constraint Ma > (Ma)w sets an extremely stringent upper bound on m 
and thus forbids the wide (empty) region on the right, ii) Si = 0,S 2 = —0.2 (Fig. 13a): here 
one still has K 2 < 0, J 2 > 0, but \K 2 \ is smaller than in the previous case, so the constraint 
Ma > {Ma)w is still very effective but less compelling than in the case (i). Also, the role of 
the m& and the b — > 57 conditions is more significant here, iii) Si = 0.7,52 — 0.4 (Fig. 14a): 
here K 2 > 0, J 2 < 0, Ma > (M^)^ gives a lower bound on mi/ 2 and the /1 2 > /i 2 b condition 
provides the frontier of the empty domain on the right. 

The (7711/2, /i) representations for large tan /3 and for the representative Si points discussed 
above are displayed in Figs.l2b-14b. We start from the universal case of Fig. 12b. Here we 
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expect gaugino-dominated configurations. However, because the values of m are strongly 
limited from above (see Fig. 12a), we have the extremely correlated states shown in Fig. 12b. 
In the case of Fig. 13b one has J 2 > 0, and gaugino-dominated states occur. No strong 
(mi/2, /i) correlation shows up in this case. The opposite case, J 2 < 0, is shown in Fig. 14b, 
where higgsino-dominated configurations appear. 



It is worth adding a few comments about the examples of Figs. 11 and 12, where the 
physical regions in the (m^, /i) plane show a very pronounced correlation in the two 
variables. This feature occurs whenever l^l^o ^ ^i m i/2' *- e -' whenever m is severely 
bounded from above and/or |J 2 | is very close to zero. As far as the values of {Jzl/Ji are 
concerned, we notice that in the universal case (see Fig.2c), except for small values of tan j3, 
J 2 /Ji ~ 0.04 (in fact, for tan/5 > 4, J x ~ -a 2 ~ 2.5, J 2 ~ -b 2 ~ 0.1). Thus for $ = a 
strong (m.1/2, //) correlation occurs whenever tuq < 0(m!/ 2 ). This happens in the example 
of Fig. 12, where m is severely bounded by the M4 > {Mj^ib condition, and in the case of 
Fig. 11, where the correlation is enforced by a very small value of J 2 - J2 — 0.06. A (mi/2, //) 
correlation is also exhibited in Fig. 9b for the range mi/2 > 150 GeV, where mo is bounded 
by the cosmological constraint. 

In general, we do not consider these physical regions with a strong (muz, fi) correlation 
as unnatural, since they are usually realized without much tuning. We recall that the size 
of the coefficients J\ and J 2 is dictated by the RGE's with their intrinsic cancellations, 
and that one naturally has J\ = 0(a few), J 2 = O(0.1 — 0.01). As we have seen, these 
properties, combined with severe upper bounds on tuq, are sufficient to generate the (rn.1/2, 
/i) correlation. 

We turn now to the A 7^ case. First we recall that A is constrained in the range 
\Aq\ < 3 from the absence of charge and color breaking [[34]. Thus, allowing A 7^ does 



not change essentially the general picture previously discussed. The previous scenarios still 
occur, but at different points in the parameter space. Two specific comments are in order 
here: i) independently of its sign, A disfavours the changeover from gaugino dominance 
to higgsino dominance in the neutralino composition, ii) a negative A reduces the value 
of Ma as compared to the A = case, and so either provides a light A boson (and hence 
interesting phenomenology) or enforces a more stringent constraint on the parameter space. 



VIII. NEUTRAL HIGGS BOSONS 



Neutralino direct detection, to be discussed in the next Section, is based on neutralino- 
nucleus scattering. In this process, exchanges of neutral Higgs bosons play a dominant role, 
provided the Higgs masses are not too heavy. It is convenient to recall here some relevant 
properties of the couplings of \ with matter via Higgs exchange. As was already mentioned 
in the Introduction, the two Higgs isodoublets H 1 , H 2 yield 3 neutral Higgs mass eigenstates: 
one CP-odd (A) state, whose mass Ma is given by expression ([□]) and two CP-even states 
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(of masses \ //,. Mh, Mh < Mh), which are obtained from H®, by a rotation through an 
angle a 

H = cos a H 1 + sin a if ° 

h = - sinaH^ + cosaH^ . (18) 



It is important to notice here that a depends very sensitively on Ma, being very close to 
zero for tan/3 > 4 and rising very fast to tt/2 for Ma < OiMz) (see Fig. 15). 



The angle a plays a crucial role in determining the size of the neutral h, if-quark cou- 
plings. Here, as we are interested in ^-nucleus scattering, we discuss explicitly only the 
couplings involving the CP-even states, since h, H are dominant compared to A. The low- 
energy neutralino-quark effective Lagrangian generated by Higgs exchange may be written 



as follows 35 



£ e fr = V2G F — 5^- 

m h,H 



(19) 



Here F^h is the ratio of the Higgs-neutralino coupling to the 5*^7 (2) gauge coupling, which 
depends on the composition of x 

Fh = 0,2 (et3 sin a + 0,4 cos a) 

Fh = 0.2(0.3 cos a — CI4 sin a) (20) 



and the k q are given, for the up-type quarks and the down-type quarks respectively, by 

H h 

k u sin a I sin f3 cos a/ sin (5 

kd cos a I cos (3 — sin a/ cos (3. (21) 



Note that, in general, since tan (3 > 1, the strength of the coupling to the down-type quarks 
is bigger than the one to the up-type quarks, and £ e ff usually gets a sizeable contribution 
when the h boson is exchanged (h is lighter than H and is therefore favored because of the 
propagator denominator in Eq. ([l9|) ) and when a ~ tt/2, i.e., when Ma < 0(Mz)- When 
this regime does not apply, the size of C e s is much suppressed. 

The cross section for elastic neutralino-nucleus scattering which follows from the effective 
Lagrangian ( |T9"[ ) will be given in Sect. IX. B. 



15 



IX. DIRECT DETECTION 



Much experimental activity is under way in the direct search for neutralino dark matter 
and the perspectives for significant improvements in experimental sensitivities are encourag- 
ing [|36|| . In this class of experiments, a relic neutralino would be detected by the amount of 
energy released by its elastic scattering off nuclei in an appropriate apparatus. A signature 
would be provided by a yearly modulation of the signal, whose observations would require 
high statistics and extremely good stability in the detector response. Here we evaluate the 
event rates for this process extending previous analyses to the non-universal 5 j 7^ case. 
Various materials are being used in the current experiments and others are under investi- 
gation for future detectors. In this paper we analyse two of the most interesting materials: 
Ge (in its natural composition) []3"T|-pD|] and 129 Xe [HO]. 



A. Differential rates 



The nuclear recoil spectrum may be evaluated from the expression 

^ = E^%^)/(^) (22) 



where 



Ro, i = N T ^a i <v> . (23) 
m x 



In Eqs.(f2lq)-(|23|) we use the following notations: the subscript i refers to the two cases of 
coherent and spin-dependent effective interactions, is the number of the target nuclei per 
unit of mass, p x is the local neutralino matter density, and E R is the nuclear recoil energy 
given by E R = m 2 ed v 2 {l — cos 6*) /m^, where 9* is the scattering angle in the neutralino- 
nucleus center-of-mass frame, is the nuclear mass, m re d is the neutralino-nucleus re- 
duced mass and v is the relative velocity. The maximum value of E R is E^ ax = 2m 2 ed v 2 jm^. 
Returning to (^2|-[23|), F(E R ) denotes the nuclear form factor, and Oi is the (coherent /spin- 
dependent) neutralino-nucleus cross section. The factor I{Er) is given by 

< V 2 > /""max f( y) 

I(E R ) = ^— / dv^ (24) 

< v > (E R ) V 



where f(v) is the velocity distribution of neutralinos in the Galaxy, as measured in the 
Earth's rest frame, and v min (E R ) is given by v min (E R ) = (m N E R / (2m 2 ed )) 1 ^ 2 . The averages 
appearing in Eqs.(|22])-(|2"4|) denote averages over the velocity distribution in the Earth's rest 
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frame. An explicit formula for I(Er) in the case of a Maxwellian velocity distribution may 
be found in Ref . f| . 

The differential rates to be discussed below will be expressed in terms of the electron- 
equivalent energy E ee rather than in terms of Er. These two variables are proportional: 
E ee = QEr where Q is called the quenching factor: typical values of Q will be discussed 
shortly. 



B. Neutralino nucleus elastic cross sections 



The total cross sections for neutralino-nucleus elastic scattering have been evaluated fol- 
lowing standard procedures PJl], |35] , [41| ,fl2]j . Here we only summarize some of the main prop- 



erties. Neutralino-quark scattering is described by amplitudes with Higgs-boson exchanges 
and Z-boson exchange in the t-channel, and by amplitudes with squark exchanges in the 
s- and u-channels. The neutral Higgs bosons considered here are the two CP-even bosons: 
h,H and the CP-odd one: A, whose couplings were previously discussed in Sect. VIII. 

The relevant properties for these amplitudes are: 1) Higgs-boson exchanges contribute 
a coherent cross section which vanishes only when there is no zino-higgsino mixture in the 
neutralino composition [|35[] , 2) Z-boson exchange provides a spin-dependent cross section 



which receives contributions only from the higgsino components of x> 3) squark exchanges 
contribute a coherent cross section (due to zino-higgsino mixing) as well as a spin-dependent 
cross section (due mainly to the gaugino components of x) ||41|| . As examples we recall here 
only the expressions for the coherent cross section due to the exchange of a Higgs boson (h 
or H) and the spin-dependent one due to Z exchange. 

The former cross section is easily evaluated from the effective Lagrangian of Eq. flT9| ) |35j 

8G 2 m 2 

a cu = — ^ai iff my 2 (25) 

71 m h,H 

where A is the nuclear mass number and ol^h is given by 

a h , H = F h , H I , I = Y,k q m q (N\qq\N). (26) 



The quantity / may be expressed conveniently in terms of the irN sigma-term a w N arid of 
a parameter a which is related to the strange-quark content of the nucleon y by 

/ // w < N\ss\N > 

a = y{m s (m u + m d )) , y = 2———— ~~T7Tjrr~ ■ 27 

< N\uu + dd\N > 

One has 
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k u g u + k d g d 



(28) 



where 
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9u = 2Y ( m ^ + — - ao-TT^ 



23 



25 



9d= 27 \ N + ~4 <J7tN + Y a<JnN 



(29) 



Unfortunately, the values of both the quantities y and a n N are somewhat uncertain. Here, 
for y we use the central value of the most recent evaluation: y = 0.33 ±0.09, obtained from a 
lattice calculation [43|. For cr^N, which is derived by phase-shift analysis and dispersion re- 
lation techniques from low-energy pion-nucleon scattering cross-sections [0,0 , we employ 
the value of Ref. gj|: cr wN = 45 MeV. We then find the results: g u = 123 MeV, g d = 288 
MeV (we use 2{m s /{m u + m d )) = 29 [f5|). We note that these values further reinforce the 
role of the down-type quarks as compared to the up-type ones. 



We point out that the Higgs-nucleon couplings for nucleons bound in a nucleus may be 
renormalized by the nuclear medium. As a consequence, the strength of I might in principle 



be reduced to some extent [47|. However, this effect is neglected here. 



Now let us turn to the spin-dependent cross section due to Z exchange. This may be 
cast into the usual form |41|J42] 



0"SD 



TT 



a 3 — a 



2\2 2 , 
4 J m redS 



Y^T 3L , q Aq) 2 \ 2 J(J+l) 



(30) 



In this paper we use this formula for 73 Ge (this isotope is present at the level of 7.8 % in 
the natural composition of Ge) and to 129 Xe. For these nuclei we employ the values of A 2 
obtained in the odd-group model 0], where only the odd nuclear species in odd-even nuclei 



are explicitly taken into account. The Ag's in Eq.(30) denote the fractions of the nucleon 
spin carried by the quarks q in the nucleon of the odd species, and the T 3 l i(? 's stand for the 
third components of the quark weak isospin. The values for the Ag's are taken from Ref. 



It is worth noticing that the event rates for neutralino direct detection with the materials 
considered here are largely dominated by coherent effects in most regions of the parame- 
ter space. In the small domains where spin-dependent effects dominate over the coherent 
ones the total rates are usually too small to allow detection. The experimental strategy of 
employing materials enriched in heavy isotopes of high spin is interesting for a search for hy- 
pothetical dark matter particles which interact with matter via substantial spin-dependent 
interactions. However, this approach does not appear to be very fruitful for neutralinos. 

One more ingredient which enters the event rate in Eq.(|22|) is the nuclear form fac- 
tor, which depends sensitively on the nature of the effective interaction involved in the 
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TABLE I. Characteristics of some current experiments. In the second column is reported 
the quenching factor Q, in the third column the electron-equivalent energy at threshold, in the 
fourth the square of the form factor at threshold, and in the last column the present experimental 
sensitivity. 



Nucleus 


Q 


E%(KeV) 


F 2 (E% ) 


evts/(Kgd KeV) 


Ge @ 


0.25 


2 


0.87 


3.0 


Ge H 


0.25 


12 


0.41 


0.2 


Xe @ 


0.80 


40 


0.07 


0.8 



neutralino-nucleus scattering. For the coherent case, we simply employ the standard pa- 
rameterization [|9 



F(E R ) 



qr 



• i 2 2 

as q 



(31) 



where q 2 =| q \ 2 = 2m^E R is the squared three-momentum transfer, s ~ 1 fm is the thickness 
parameter for the nucleus surface, ro = (r 2 — 5s 2 ) 1 / 2 , r = 1.2 A 1 / 3 fm and ji(qro) is the 
spherical Bessel function of index 1. 

The form factor in Eq.([H|) introduces a substantial suppression in the recoil spectrum 
unless qr < 1. A noticeable reduction in dR/dE R may already occur at threshold E R = 



E R = E^/Q when rQ\j2mNE R h is not small compared to unity. The actual occurrence 
of this feature depends on parameters of the detector material: nuclear radius, quenching 
factor, threshold energy E^. The values of these parameters for the nuclei considered in 
this paper are reported in Table I |36|j3^ -^0| , and the values of F 2 (E R ^ 1 ) calculated from 



Eq.fplD are given in this same Table. Since we consider in this paper mainly the value of 
the differential rate near threshold, F 2 (E R ^ 1 ) is the most relevant quantity. We see from the 
values in the Table that the reduction introduced by the form factor is moderate in Ge, but 
quite substantial in 129 Xe. 

In general, for the spin-dependent case there are no analytic expressions for the form 
factors. However, numerical analyses have been performed for a number of nuclei. The 
general feature is that these form factors have a much milder dependence on E R as compared 
to the coherent ones, because only a few nucleons participate in the neutralino-nucleus 
scattering in this case. In our evaluations we use the results of Refs. for 131 Xe and 

73 Ge respectively. 
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C. Local Neutralino Density 



cm 



We denote the local halo density by pi, for which we use the estimate pi = 0.5 GeV 
~ 3 fn| . For the value of the local neutralino density p x to be used in the rate of Eg. (p3|) . 



for each point of the model parameter space we take into account the relevant value of the 
cosmological neutralino relic density When Vt x h 2 is larger than a minimal (Qh 2 ) m m required 
by observational data and by large scale structure calculations we simply put p x = pi. When 
Q x h 2 turns out less than (Qh 2 ) min , the neutralino may only provide a fractional contribution 
fl x h 2 /(flh 2 ) min = £ to Qh 2 ; in this case we take p x = pi£. The value to be assigned to 
(Qh 2 ) min is somewhat arbitrary. Here we set it equal to 0.1. 



It is worth remarking here that, due to this scaling procedure, for the direct detection 
rate one has: i) i?o,i oc piOi for Q x h 2 > (Qh 2 ) m i n and ii) _Ro,i oc pi^di oc piOij < cr ann t> >; nt for 
Q x h 2 < (Qh 2 ) min . Thus the rate i? ,j is large in the regions of the parameter space where <7j 
is large. This is trivial in case i), but it is also true in case ii), since when <jj is large also 
<r ann increases but in such a way that usually the ratio o"j/<T ann increases too. Because of the 
relation VL x h 2 oc< cr. Ann v > in J it follows that Rq^ is large for neutralino configurations with 
modest values of the relic abundance, and vice versa. 



D. Results for detection rates 



The most significant quantity in comparing experimental data and theoretical evaluations 
for direct detection is the differential rate dR/dE ee = (dR/dEji) /Q (with dRjdE^ defined 
in Eq.(^)) rather than the total rates, obtained by integration over wide ranges of E ee . 
By using the differential rate instead of the integrated ones, one obtains the best signal- 
to-background ratio. Note that the experimental spectra, apart from an energy interval 
around threshold, usually show a very flat behaviour, whereas signals for light neutralinos 
are decreasing functions of the nuclear recoil energy. 

A complete procedure would then be to compare the experimental and theoretical rates 
over the whole E ee range. However, to simplify the presentation here, we give our results 
in terms of the rate integrated over a narrow range of 1 KeV at a specific value of E ee , the 
one which appears the most appropriate for each experiment: typically it corresponds to a 
point close to the experimental threshold. To be definite we consider the following cases: 



i) Ge (natural composition). Among the various running experiments [3b|, we select 
the two which, at present, appear to provide the most stringent limits: a) Caltech-PSI- 

3 events/ (Kg day KeV); b) Heidelberg- 
0.2 events/(Kg day KeV). Correspond- 
ingly, for Ge we have evaluated our rate by integrating dR/dE ee over the range (2-3) KeV 
for experiment a) and over (12-13) KeV for experiment b). It turns out that the case b) 



Neuchatel @ with E. 
Moscow |39ll with E. 



th 

ee 



tli 



= 2 KeV, differential rate 
12 KeV, differential rate 
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provides the most stringent bound also for light neutralinos. 



ii) Xe. In this case, taking into account the features of the DAMA experiment |40 



we 



have considered the rate R integrated over the range (40-41) KeV. 



Our results are shown in Figs. 16-19. Figs. 16-18 report the rate for a Ge detector for the 
regions of the parameter space which are depicted in Figs. 12-14, respectively. In parts (a) 
and (b) of each figure, R is displayed in the form of a scatter plot, in terms of m x and of 
the relic abundance, respectively. The horizontal line denotes the present level of sensitivity 
in the Heidelberg-Moscow experiment. We notice that, in all cases shown in these figures, 
the experimental sensitivity is already, for some configurations, at the level of the predicted 
rate. Some points of the supersymmetric parameter space, denoted by filled squares in 
Figs. 12-14, are even already excluded by present data. The exploration potential of this 
class of experiments as the sensitivity is improved is apparent from these figures. Fig. 19 
shows the rate R for 129 Xe for the region of the parameter space displayed in Fig. 13: again 
the horizontal line gives the present experimental sensitivity. A comparison of Fig. 19 with 
Fig. 17 shows that the Ge experiments are currently more effective. However, it has to be 
noticed that experiments with liquid Xe may become extremely competitive in the future 



A few more remarks are in order here: 



i) The cases displayed in Figs. 16-18 present the common feature of providing fair chances 
for direct detection. This is not a surprise, since these representative points all belong to the 
category of configurations with small values of M4. As was stressed before, once we move 
away from these appealing physical regions of the neutralino parameter space, the rates for 
direct detection may fall far below (by many orders of magnitude) the detection sensitivities 
(present or future). This unfortunate situation occurs, for instance, typically as we move 
towards smaller values of tan /3. However, one should keep in mind that the regime of very 
large tan f3, where signals may be sizeable, represents a very interesting scenario, deserving 
much attention and exploration. In fact this is one of the two options, very small or very 
large tan/3, which seem to fit low-energy phenomenology at the best [ p2]| . 

ii) The scatter plots in parts b) of Figs. 16-19 show explicitly a property previously men- 
tioned in Sect.IX.C, namely that the scaling procedure adopted to evaluate the neutralino 
local density implies a K-Q x h 2 correlation. Configurations which provide a measurable R 
usually entail a low Q and viceversa. Only in a few cases the neutralino may be detectable 
by direct detection and also provide a sizeable contribution to Q. 
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X. CONCLUSIONS 



In the present paper we have discussed some possible scenarios for neutralino dark matter 
which originate from the relaxation of the assumption of strict universality for soft scalar 
masses at M GUT . 

This approach derives from the general consideration that many crucial theoretical points 
entering not only grand unified and supersymmetric theories, not to mention the Standard 
Model, are far from being understood and/or verified. For this reason, any new theoretical 
assumption has to be fully scrutinized. This is even more important because new assump- 
tions in supersymmetric models are often introduced not because of solid arguments, but 
rather for the sake of simplicity and for the need to reduce the large number of free param- 
eters that would otherwise prevent any firm prediction. 

In our work we have discussed different scenarios, by considering various physical con- 
straints in a sort of hierarchical order, giving top priority to the requirement of radiative 
EWSB, implemented with a no-fine-tuning criterion, and to the cosmological relic neu- 
tralino density constraint. Some other assumptions, often introduced in the literature, have 
been relaxed in our work. This is in particular the case for universality in the soft scalar 
masses. However, it has to be remarked that the type of departure from universality that 
we have considered in our paper is far from being the most general one, as was noticed in 
Sect. II. In particular, it only refers to the Higgs masses, and not to the sfermion masses. 

The implications of the various scenarios on neutralino relic abundances and rates for 
detection rates have been analysed, and the impact of a non-universality in has been dis- 
cussed for the whole range of tan/3. We have shown that the departure from m universality 
is particularly interesting in two respects: 

i) Small values of Ma are allowed: this has in itself the dramatic consequence for direct 
detection of generating a large value for the angle a and large couplings to matter of the 
lightest neutralino x- 

ii) Higgsino or mixed higgsino-gaugino configurations appear for all tan/3: this contrasts 
with the pure gaugino configurations favoured by strict m universality. 
Consequences of such a departure from universality on the size of the neutralino relic abun- 
dance have been analysed for both large and small values of tan/3. It has been shown that, 
because of the previous properties, deviations from universality may reduce the value of 
Q x h 2 . 

The predicted rates for direct detection has been analysed in detail and compared with 
current and foreseen experimental sensitivities. The role of the previous properties in opening 
interesting perspectives for this kind of search has been elucidated. We find that presently- 
running experiments are already impacting interesting regions of the neutralino parameters 
space in some of the non-universal scenarios discussed here. 
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Figure Captions 



Figure 1 - The (mi/ 2 , //) plane for tan/? = 8. The lines of constant m x = 30 GeV, 60 
GeV, 90 GeV are displayed as dashed lines. The lines of constant P = 0.1, 0.5, 0.9, 0.99 are 
shown as solid lines. The dotted region denotes the domain forbidden by present LEP data. 



Figure 2 - Coefficients of the polynomial expressions (9-11) as functions of tan/3: 

a) coefficients of Mj^: a\ solid line, b\ dashed line, c\ dot-dashed line and d\ dotted line, 

b) coefficients of Mjj 2 : a 2 solid line, b 2 dashed line, C2 dot-dashed line and d 2 dotted line, 

c) coefficients of fi 2 : J\ solid line, J2 dashed line, J3 dot-dashed line and J4 dotted line, 

d) coefficients of M\: K\ solid line, K 2 dashed line, K 3 dot-dashed line and K4 dotted line. 



Figure 3 - In the (S 2 , Si) plane, the lines where J2 = at fixed tan/5 are displayed: 
tan/? = 53, solid line; tan/3 = 8, dashed line; tan/3 = 3, dot-dashed line; tan/3 = 2, dotted 
line. 



Figure 4 - In the (S 2 , Si) plane, the lines where K 2 = at fixed tan/3 are displayed: 
tan/3 = 53, solid line; tan/3 = 40, dashed line; tan/3 = 8, dot-dashed line. 



Figure 5 - Graph of Ma as a function of tan/3 at the representative point m = 50 
GeV, m V 2 = 200 GeV. 



Figure 6 - Scatter plot representing VL x h 2 as a function of m x for tan/3 = 8, Si — and 
S 2 = 0. Parameters are varied on a linear equally-spaced grid over the ranges: 10 GeV < 
m < 2 TeV, 45 GeV < m 1/2 < 500 GeV. 



Figure 7 - The same as in Figure 6, but with Si = —0.2 and S 2 = 0.4. 
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Figure 8 - Scatter plot representing Q x h 2 as a function of m x for tan/3 = 53, 8± = 
and 62 = 0. Parameters ranges are as in Figure 6. 



Figure 9 - a) Parameter space in (m 1( / 2 , m ) plane for tan/? = 8, Si = and 5 2 — 0. 
Empty regions are excluded by: i) accelerator constraints, ii) radiative EWSB conditions, 
iii) neutralino is not the LSP. Dots represent the region where Q x h 2 > 1. Regions with 
crosses are excluded by b — > 57 and m& constraints. In the regions denoted by squares, 
0.01 < Q x h 2 < 1. The region without fine-tuning is inside the box bounded by dashed 
lines. 

b) Parameter space represented in the (myz, £t) plane. Solid lines correspond to the extreme 
values of mo- Notations are the same as in a), but crosses are omitted here. 



Figure 10 - a) The same as in Figure 9, but with Si = —0.2 and S2 = 0.4. In the regions 
denoted by diamonds fl x h 2 < 0.01. 

b) Notations are the same as in a). Note that the domain where the neutralino is the dark 
matter particle with 0.01 < Q x h 2 < 1 has shifted to the higgsino-dominated region. Crosses 
are omitted here. 



Figure 11 - The same as in Figure 9, but with 8\ = —0.8 and 5 2 = 0.2. 



Figure 12 - a) Parameter space in the (m^, m ) plane for tan/3 = 53, Si = and 
S 2 = 0. Filled squares denote configurations excluded by direct detection with a Ge detector 
p9| . Other notations are as in Figure 9. 

b) Parameter space represented in the (m 1( / 2 , fi) plane. Notations are the same as in a). 



Figure 13 - The same as in Figure 12, but with Si — and S2 = —0.2. 



Figure 14 - The same as in Figure 12, but with Si = 0.7 and S2 = 0.4. 



Figure 15 - a) Graphs of Ma and as functions of tan/3, 
b) Mixing angle a as a function of tan /3. 
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Here the representative point is m = 100 GeV, = 200 GeV. 



Figure 16 - Scatter plot of the rate for direct detection with a Ge [B5] detector for 
tan/3 = 53, Si = and S2 = 0, as a function of m x (a) and as a function of Q x h 2 (b). 
Parameters are varied on a linear equally-spaced grid over the ranges: 10 GeV < tuq < 
2 TeV, 45 GeV < m 1/2 < 500 GeV. 



Figure 17 - The same as in Figure 16, but with Si — and 5 2 = —0.2. 



Figure 18 - The same as in Figure 16, but with Si = 0.7 and 82 = 0.4. 



Figure 19 - Scatter plot of the rate for direct detection with a Xe 0] detector for 
tan/3 = 53, Si = and S 2 = —0.2, as a function of m x (a) and as a function of £l x h 2 (b). 
Parameters are varied on a linear equally-spaced grid over the ranges: 10 GeV < m < 
2 TeV, 45 GeV < m 1/2 < 500 GeV. 
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